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ABSTRACT 



Aims. The aim of this paper is to demonstrate the effect of turbulent background density fluctuations on flare-accelerated electron 
transport in the solar corona. 

Methods. Using the quasi-linear approximation, we numerically simulated the propagation of a beam of accelerated electrons from 
the solar corona to the chromosphere, including the self-consistent response of the inhomogeneous background plasma in the form 
of Langmuir waves. We calculated the X-ray spectrum from these simulations using the bremsstrahlung cross-section and fitted the 
footpoint spectrum using the collisional "thick-target" model, a standard approach adopted in observational studies. 
Results. We find that the interaction of the Langmuir waves with the background electron density gradient shifts the waves to a higher 
phase velocity where they then resonate with higher velocity electrons. The consequence is that some of the electrons are shifted 
to higher energies, producing more high-energy X-rays than expected if the density inhomogeneity is not considered. We find that 
the level of energy gain is strongly dependent on the initial electron beam density at higher energy and the magnitude of the density 
gradient in the background plasma. The most significant gains are for steep (soft) spectra that initially had few electrons at higher 
energies. If the X-ray spectrum of the simulated footpoint emission are fitted with the standard "thick-target" model (as is routinely 
done with RHESSI observations) some simulation scenarios produce more than an order-of-magnitude overestimate of the number of 
electrons > 50keV in the source coronal distribution. 

Key words. SumCorona - Sun:Flares - Sun: X-rays, gamma rays 



1. Introduction 

The unprecedented RHESSI observations of solar flare hard 
X-rays (HXRs, typically > 20keV) has forced us to consider 
mechanisms in addition to the traditional collisional view of 
coronal electron transport. This standard approach is of an as- 
sumed power-law of electrons above a low-energy cut-off that 
propagate downwards, losing energy through Coulomb collision 
with the "cold" background plasma (whose energy is consid- 
erably lower than that of the electrons in the beam). The elec- 
trons will eventually stop once they have reached the higher 
density chromosphere (the "thick-target"), emitting X-rays via 
bremsstrahlun g and heating the plasma. The "cold thick-ta rget" 
model CTTM (Brownl ll97lHSvrovatskii & Shmelevalll972l) has 
proved popular because it provides a straightforward relation- 
ship between the bright X-ray emission from the chromospheric 
footpoints and the source coronal electron distribution. However, 
many RHESSI observations a re not consistent with the CTTM 
or pro duce challenging results (Holman et al. 20lTl lKontar et all 
I201 ll) . which demonstrates that essential physics is missing from 
the standard model. 

One aspect missing from the CTTM are non-collisional pro- 
cesses such as wave-particle interactions. The self-consistent 
generation of Langmuir waves by the electron beam is one 
such process that is thought to be the related to the decimet- 
ric rad io emission in reverse-slope Type III bursts in some 
flares (iTarnstrom & Zehntnerl Il975t lAschwanden et all Il995t 



iKlein et alj[l997l; lAschwanden & Benzlll997h . We have previ- 
ously shown that including Langmuir waves helps in allevi- 
ating the discrepancies between the CTTM and RHESSI ob- 
servations. For instance, the CTTM predicts a "dip" to ap- 
pear in the flare electron spectrum between the thermal com- 
ponent and just before the turnover in the electron beam spec- 
trum. Observationally this has not been confirmed and we 
showed that the growth of Langmuir waves flattens the elec- 
tron spectrum at lower energies, maintaining a negative gra- 
dient between the therm al and non-thermal spectral compo- 
nent dHannah et al.l |200"9|) . Another observational challenge is 
the difference in the HXR spectral indices found between the 
footpoint and coronal sources, which the CTTM predicts to be 
Ay = 2, yet RHESSI's imagi ng spectroscopy of some flares has 
found values with Ay > 2 (Emslie et al. 2003; Battagl ia & Benz 
2007; ISaint-Hilaire et al.ll2008HBattaglia & BenzH2008tlSu eTaL 
2009). We found that the flattening of the electron distribution 
as it propagates from the corona to the chromosphere due to the 
generation of Langmuir waves can p roduce the observed larger 
differences (Hannah & Konta 
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Wave-particle interactions can produce some observational 
features of flares better than the CTTM but the flattening of the 
electron distribution to lower energies through Langmuir wave 
growth produces far fainter HXR emission. This means that a 
higher number of electrons need to be accelerated in the corona 
for the simulations including wave-particle interactions to pro- 
duce a similar magnitude of HXRs to the standard collisional 
approach. Compounding this problem further is that the CTTM 
itself needs a large number of electrons to be accelerated in the 
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corona, which conflicts with the maximum resupply rate of elec- 
trons in the coronal acceleration region in some models. 

The Langmuir waves themselves might provide a solution to 
this problem as they can be scattered or refracted when interact - 
ing with an inhomogeneous background pla sma dRvutovll 19691). 
which can result in el ec tron acceleration (e.g.lMelrose & Cramer 
ll989tlKontdl2001allrllR~eid & Kontad 12010) and increased X- 
ray emission (iKontar et al.l2012h . The core idea is that the waves 
can be shifted to a lower wavenumber (higher phase velocity) by 
interacting with the density gradient in the background plasma, 
which then resonates with electrons at higher velocity. Although 
this can happen in the opposite direction (with the waves shifted 
to higher wavenumber), the falling power-law electron distribu- 
tion always means that this effect has the strongest consequences 
for the re-acceleration of electrons to higher energies. Recently, 
the role of the non-uniform plasma has been studied fo r an in- 
terplanetary electron beam dReid & Kontarll2010l 1201 2l) and in 
the statio nary (no spatial e volution) case for solar flares in the 
corona (IKontar et al.l2012l) . It was found that this can lead to ad- 
ditional electron acceleration. In these studies different forms of 
the inhomogeneity in the background plasma were considered, 
including a Kolmogorov-type power-density spectrum of fluctu- 
ations, which imitates the spectra expected from low-frequency 
MHD-turbulence. 

In this paper, we demonstrate the consequences of this self- 
consistent treatment of electron beam-driven Langmuir waves 
propagating through an inhomogeneous background plasma. 
This is simulated using the quasi-linear weak-turbulence ap- 
proach and is detailed in §0 The resulting electron and spec- 
tral wave density distributions for a variety of forms of the input 
electron beam are shown in $3] The mean electron (deducible 
from observations) and X-ray spectra are obtained for these sim- 
ulations and the latter are fitted as if they were observations, us- 
ing the standard CTTM approach. This allows us to determine 
the discrepancy between the CTTM-derived and true properties 
of the source electron distribution in 33.11 

2. Electron beam simulation 

Following the previously adopted approach (lHannah et al.B2 009; 
lHannah &Kontatll20 111) , we simulated a ID velocity (v ~ V\\ » 
Vj.) electron beam f(v, x, t) [electrons cm 4 s] from the corona 
to the chromosphere, that self-consistently drives Langmuir 
waves (of spectral energy density W(v, x, t) [erg cm -2 ]). 
This weakly turbulent description of quasi-linear relaxation 
dVedenov & Velikhov|[T9 63: Dru mmond & Pines 19 641 iRvutov l 
1969 1 ; lHamilton & Petrosian|[l987l:lKontaril2001 at lHannah et alj 
2009) is given by 

df | v df _ An 2 e 2 8 (Wdf \ | ^_[L + V \ d l\ (1) 
dt dx ml dv \ v dv) dv \ v 2 v 3 dv t 

dW 3v 2 T dW v 2 dW (nu p 2 df L .. 

(2) 

where n the background plasma density, m e the electron mass 
and u> 2 = \nne 2 jm t is the local plasma frequency. The first 
terms on the right-hand side of Eqs. (QJ and (f2) describe 
the quasi-linear interaction, the other terms the Coulomb col- 
lisions yc F - 47re 4 n In A/m 2 and yc w = ne A n In A/(mlvj) 
with In A the Coulomb logarithm, Landau damping = 



■\Jn/8u) p (v/Vt) 3 exp (-v 2 /2v 2 ), and spontaneous wave emis- 
sion S = oAmv In (v/vt)/(47t«). The simulations here fea- 
ture two changes over the previous work (lHannah et al.l 12009: 
lHannah & Kontaill201 ll) . The first minor change is the inclusion 
of the diffusion in velocity space due to collisions (final term in 
the brackets at the end of Eq. (JTJ). Previously only the drag term 
was used, which described a "cold" target situation in which the 
energy in the beam electrons is considerably higher than that 
of the background thermal distribution. Including the diffusion 
term allows a more realistic treatment of electrons at energies 
closer to the thermal background, i.e. a "warm" target. The sec- 
ond, and more substantial, change is the inclusion of a turbulent 
background plasma and the effect of this density gradient on the 
plasma waves (third term on the left-hand side of Eq. ©. This 
is done using the characteristic scale of the plasma inhomogene- 
ity L' 1 _ {d(x> v ldx) <x>~ 1 — (dn/dx) ( 2n)~ l as used previously by 
IKontar (2001a) and Reid_&jConto| «HqT3). 

The electron distribution is simulated in the velocity domain 
from the 1MK background plasma thermal velocity vt up to 
115vj This range extends below the observational capacity of 
RHESSI (down to about 86eV instead of RHESSI's 3keV limit), 
well into the range where the thermal emission will dominate. 
Although the treatment of the thermal distribution is beyond 
the scope of this work, the energy range is included to demon- 
strate the possible effect these lower energy electrons have on 
the higher energy population. 

The initial electron distribution is Gaussian of width d = 2 X 
10 8 cm and a broken power-law in velocity, which is flat below 
the break. This break is effectively the low-energy cut-off. We 
used vc ~ 9vt (Eq = 7keV), and the power-law above it has an 
index of 2<5b (hence a spectral index of <5b in energy space), i.e. 

/(v,x, f = 0)oc nb exp^-j| (v/vcrMb Vc < v<V0; O) 

where vo = 90vt is the maximum initial beam velocity, «b = 
f /dv is the electron beam density above the break/low-energy 
cut-off. For the simulations presented in this paper we used a 
beam density above the break of «b = 10 8 crrT 3 . At the start 
of the simulations this beam was instantaneously injected at a 
height of 40Mm above the photosphere and was not replenished, 
with the spatial grid extending f rom 52Mm d own to 0.3Mm. 

A finite difference method riKontaril200l"cl) is used to solve 
Eqs. (QJ and ©, and the code is modular which makes it easy to 
consider the effects of the different processes. We consider three 
distinct simulation setups within this paper, namely: 

- B: beam-only: We only consider the propagation of the 
electron distribution subject to Coulomb collisions with the 
background plasma, similar to the CTTM. Specifically, we 
only solve Eq. ([TJ, ignoring the quasi-linear term (first term 
on the right-hand side). 

- BW: Beam and waves: We consider the propagation of 
the electron distribution and the self-consistent driving of 
Langmuir waves but without the plasma inhomogeneity 
term. This is solving both Eqs. ([TJ and (f2]i but without the 
last term on the left-hand side of Eq. (0, i.e. ignoring L _1 . 

- BWI: Beam, waves and inhomogeneity: We consider the 
propagation of the electron distribution, the self-consistent 
driving of Langmuir waves, and the wave interaction with the 
inhomogeneous background plasma. This is the full solution 
to both Eqs. (QJ and ©. 

These three different simulation setups allow us to investigate 
electron transport due to collisions (B, akin to the CTTM), colli- 
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Fig. 1. Background plasma density profile n. The inset view is of 
a zoomed portion of the density profile, showing that the density 
fluctuations (1000 of them with wavelengths between 10 3 < A, < 
10 6 cm and amplitude of 1 %) have been added to the background 
profile. 
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Fig. 2. Magnitude of the characteristic scale of the plasma inho- 
mogeneity Zr 1 as a function of height above the photosphere. It 
is shown for two different wavelength ranges, 10 3 < Aj < 10 6 
cm and 10 3 < Aj < 10 8 cm, using amplitudes of 1% and 10%. 



where Aj are the wavelengths of the density fluctuations, and C 
is a normalisation constant used to control their amplitude via 
An)/( n) = (C 2 £f //2) I/2 , as implemented in iReid & Kontail 
2010). Here N — 1000 fluctuations are added to the background 
density profile with random phases </>,- between and 2n and 
wavelengths of either 10 3 < A t < 10 6 cm and 10 3 < A t < 
10 8 cm chosen randomly in logarithmic space. We investigated 
two wavelength ranges. One extends to the Mm range, though in 
both cases it is smaller than the whole simulation region, choos- 
ing C to achieve amplitudes of either 1% or 10%. The inset plot 
in FigureQ]shows the fluctuations for the 10 3 < At < 10 6 cm and 
1% case, which are also present in the main plot. The most im- 
portant aspect of the fluctuations is not the wavelength range or 
amplitude used, but the resulting magnitude of the density gradi- 
ent, which influences the waves via the characteristic scale of the 
plasma inhomogeneity L = 2n(dn / 'dxY 1 in Eq. d2]i. To calculate 
this the density gradient of the fluctuations has to be analytically 
found so that they are accurately included, i.e. 



^ = ^£ + C V ?p-A m sinfc^/) 
dx dx j~t \ A i / 



+2nonA^ 1 1 cos I h > 

1 Ai 



2nx 



(5) 



The resulting L for the four different configurations of the fluc- 
tuations is shown in Figure [2] 



2.2. Simulated X-ray and mean electron spectrum 

From these simulations we can compute the X-ray spectrum 1(e) 
and the mean electron flux spectrum (nVF(E)), deducible from 
the observed X-ray spectrum. For the X-ray spectrum 1(e) we 
used 



f(v, x,t) 
n(x) Q(e, E) 



dEdx, dt 



(6) 



where A is the area of the emitting plasma (which we took 
to be the square of the full width at half-maximum of the 
Gaussian spatial distribution, i.e. 8(ln 2d 2 )), and Q(e, E) is th e 
bremsstrahlung cross-section dKoch & Motzlll959tlHauglll997h . 
We calculated the mean electron flux spec trum (nVF(E)), wh ich 
is deducible from the X-ray spectrum (e.g. lBrown et al.l2006l) by 



sions and wave-particle interactions (BW), and collisions, wave- 
particle interactions and the plasma inhomogeneities (BWI). 

2.1. Background density fluctuations 

The main component of the background density no is constant 
in the corona (10 10 ctrT 3 ) and sharply rises through the transition 
region and chromosphere ( below 3Mm). It is s hown in Figure Q] 
and was used previously in lHannah & Kontarl d201 lh . The addi- 
tional components presented in this paper are the density fluc- 
tuations, which are drawn from a turbulent Kolmogorov-type 
/? = 5/3 power density spectrum, i.e. 



n — no 



1 +C |>f sin(^ + , 



(4) 



(nVF(E)) = 



4ttR 2 



EE 



n(x) 



f(v,x,t) 



dxdt, 



(7) 



where F(E)dE = vf(v)dv is the electron flux spectrum as a func- 
tion of energy, not velocity. Most of the spectra shown in l^^ 6 
summed over the whole simulation, but for the X-ray footpoint 
spectrum ( 93.11 ) the summation is over 0.3 < x < 3Mm. 

3. Simulation results 

Results from one configuration of the three simulation setups is 
shown in Figure[3] all using an initial beam of 5b — 3.5 and back- 
ground density fluctuations of 10 3 < At < 10 6 and An/n = 1%. 
Shown here are the electron f(v, x, t) and wave spectral energy 
distributions W(v, x, f), in terms of velocity-vs-distance travelled 
in each frame with time increasing from left to right. The first 
two configurations, beam-only (B, top panel) and beam and 
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Fig. 3. Evolution of the electron /(v, x, t) and wave spectral energy distributions W(v, x, t), with time increasing left to right, for 
different simulation setups (top to bottom) but all with 5b = 3.5 and rib - 10 8 cirT 3 . (Top) Electron distribution for the beam-only B 
simulation. (Middle) Electron distribution and wave spectral energy distribution for the beam and waves BW simulation. (Bottom) 
Electron distribution and wave spectral energy distribution for the beam, waves and inhomogeneity BWI simulation. 



waves (BW, middle panels), show similar results to those that 
we have previously published dHannah et al.ll2009h . Here we are 
using a higher beam density, however we have a flat (instead of 
no) electron distribution below Ec, and there are density fluc- 
tuations in the background plasma. In the beam-only case we 



see that the fastest electrons reach the lower atmosphere first, 
quickly lose energy to the high-density background plasma and 
leave the simulation grid. The bulk of the electron distribution 
takes longer to lose energy through Coulomb collisions with the 
background plasma. After 1 second in simulation time the elec- 
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tron distribution is no longer present. When the wave-particle in- 
teractions are included (BW, middle panel of Figure [3), the elec- 
tron distribution immediately flattens/widens in velocity space, 
with electrons shifted to lower energies. After t = 0.03s two 
components of Langmuir waves have developed: one at lower 
energies through the spontaneous emission, the S f term in Eq. 
(f2| and another across a wide range of velocities through the 
propagation of the fastest electrons away from the bulk of the 
distribution, the df/dx term in Eq. ©. 

The result that these simulations produce similar results to 
our previous work confirms that the density fluctuations only 
play a role once the inhomogeneity term L term is included in 
Eq. (0, which is shown in the bottom panels of Figure ([3J. The 
BWI shows a dramatic change over the other simulations, with 
streaks appearing in the wave spectral density plots because the 
waves shift to lower and higher phase velocity (or higher and 
lower wavenumber). The effect on the electron distribution is to 
pull it out in clumps across the velocity-space, which is most ev- 
ident in the t = 0.10s frame. The leading edge of the electron 
distribution is clearly pushed out to higher velocities compared 
to the other setups, i.e. at times t = 0.10,0.15s the frames of 
the electron distributions of the B (black) and BW (blue) setups 
do not extend to energies as high as the BWI (purple) setup in 
Figure [3] 

The change in energy is most evident when the spatially- 
integrated mean electron spectra (nVF(E)) are calculated for 
the simulations, as shown in Figure|4] Here all configurations of 
the simulations are shown, indicated by different coloured lines, 
and the panels show the spectral indices of the initial distribu- 
tion 5b - 2.5, 3.5,4.5, and 5.5 increasing from left to right. The 
mean electron spectrum of the simulations shown in Figure [3] 
are shown in the second plot in Figure |4] the same colours are 
used for the electron distributions in each different setup. For 
the hardest (i.e. flattest, 6b = 2.5, left panel Figure |4]i initial 
spectrum almost all different setups are substantially lower than 
for the beam-only case. Only the fluctuations with the steepest 
density gradient (10% and 10 3 < A, < 10 6 cm) produce elec- 
trons at higher energies than the purely collisional setup, but this 
only occurs at the highest energies (> lOOkeV). With steeper ini- 
tial spectral indices (larger 6\,) we find that more electrons have 
been accelerated with even lower levels of density fluctuations, 
though the inhomogeneity needs to have L _1 > 10~ 8 ctrr'. This 
may depend on the way the initial distributions were normalised. 
The same beam density above the low-energy cut-off was used 
throughout i% = J v fdv, but for steeper spectra this results in 
more electrons at energies just above Eq. Therefore there is a 
higher number (about a factor of two from Figure |4} of electrons 
with about 10 keV with the steeper spectra available to be re- 
accelerated by the shifted waves. Even for the strongest electron 
acceleration caused by the density fluctuations, it is only above 
about 20 keV that there are more electrons than in the beam-only 
setup. In all simulations where the wave-particle interactions are 
present, the Langmuir wave generation flattens the spectrum, 
which reduces the number of low-energy electrons compared to 
the purely collisional case. 

The result of the wave scattering for the HXR spectrum is 
more complicated since an electron can produce an X-ray below 
its energy. This is because higher energy electrons can travel far- 
ther into the dense regions of the lower solar atmosphere, pro- 
ducing substantially stronger HXR emission. The spatially inte- 
grated X-ray spectra for the different simulation setups and con- 
figurations are shown in Figure [5] As with the mean electron 
spectrum, the most significant changes are observed in the sim- 
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Fig. 6. Footpoint X-ray spectrum for an initial electron power- 
law index of 6b = 3.5 and the different simulation setups 
(coloured lines). Each is fitted over 20 - 100 keV using the 
f_thick2 .pro model, which is indicated by the dashed lines 
and is similar to the beam-only simulation (B, black line). 

ulations with the softest (steepest r5b = 5.5, last panel of Figure 
[5} initial electron distributions. Here the X-ray emission is up 
to several orders of magnitudes higher than the beam-only case, 
whose spectrum is flatter down to 10 keV. This trend contin- 
ues with the harder initial electron distributions producing flatter 
X-ray spectrum compared to the purely collisional case. Again 
with the hardest initial spectrum (Sb = 2.5, first panel Figure 
[5} only the strongest density fluctuations produce X-ray emis- 
sion higher than the beam-only case, but this extends to about 
30 keV, whereas in the electron spectrum it is only higher > 100 
keV. 

3.1. Fitting the footpoint X-ray spectrum 

To quantify the effect of the wave-particle interactions and den- 
sity fluctuations on the X-ray spectrum, we fitted them as if 
they were actual observations. We specifically fitted the foot- 
point X-ray spectrum, where in Eq. © we summed over the 
region 0.03 < x < 3Mm instead of the whole simulation, be- 
cause in RHESSI observations the flare spectrum is mostly dom- 
inated by the footpoint emission from the chromosphere. These 
spectra were fitted using the implementation of the CTTM in 
the OSPEX software f_thick2 .pro, an optimised v ersion of 
the rou tine by G. Holman (usage examples are given in Holman 
d2003h ). available in the SolarSoft X-ray packag^B We fitted a 
single power-law as the source spectrum, setting the low-energy 
cut-off (Ec = 7keV) and maximum energy to be the same as the 
initial electron distribution in our simulations. The fitting can be 
highly sensitive to the low-energy cut-off, so by fixing it to the 
true simulation value, we avoided this problem, as well as the is- 
sue of the missing thermal component at low energies that would 
be present in real spectral observations. We therefore have two 
free parameters in our model fit: the total number of electrons 
N(> Ec) and the spectral index r5jT of the source distribution. 
We fitted, by minimizing^ 2 , the simulated footpoint spectra over 
20 to 100 keV, the typical energy range used in RHESSI observa- 
tions, which also avoids complications of the thermal component 
at low energies and simulation edge effects at higher energies. 

1 http://hesperia.gsfc.nasa.gov/ssw/packages/xray/ 



5 




Fig. 4. Mean electron flux spectrum (spatially integrated and time averaged) for initial electron beams of power-law index 8b - 
2.5, 3.5, 4.5, and 5.5 (increasing left to right). The different colour lines indicate the different simulation setups : B (black dashed), 
BW (blue), and BWI (orange, green,purple, and red). In the last case two different wavelength ranges are used (10 3 < At < 10 6 cm 
and 10 3 < A[ < 10 s cm) with fluctuations of the amplitudes of 1% and 10%. 
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Fig. 5. X-ray spectrum (spatially integrated and time averaged) for initial electron beams of power-law index 6b - 2.5,3.5,4.5, 
and 5.5 (increasing left to right). The different colour lines indicate the different simulation setups : B (black dashed), BW (blue), 
and BWI (orange, green,purple, and red). In the last case two different wavelength ranges are used (10 3 < A; < 10 6 cm and 
10 3 < Ai < 10 s cm) with fluctuations of the amplitudes of 1% and 10%. 



The simulated footpoint spectra (colour lines) and their 
f_thick2 .pro fits (dashed lines) for the initial spectral index 
Sb = 3.5 are shown in Figure [6] for all simulation setups. In all 
cases the model fits the simulated spectra very well over the cho- 
sen energy range (indicated by the dotted vertical lines). The fit- 
ted £_thick2 .pro model should be a reasonable match to our 
simulation B and we obtain 6jj = 3.7 and N(> Ec) — 4 x 10 35 
electrons. The slight discrepancy between the fitted and true 
spectral index (3.7 vs 3.5) for the source distribution is because 
the £_thick2.pro model is steady-state and stationary where 
as our simulation includes the time and ID-spatial evolution of 
an injected (not continuous) electron beam. For the simulations 
with different initial spectral indices we again find only a small 
discrepancy to the fitted values, obtaining 2.8, 4.7, and 5.6. In all 
these B simulations we obtained the total number of electrons in 
the range of N(> E c ) = 3.5 - 4.7 x 10 35 electrons. 

We normalised all fitted results by those found for the B case. 
They are shown for the spectral indices £tt in Figure [7] and total 
number of electrons N(> Eq) in Figure|8] For the steepest initial 
distribution (5b = 5.5) all fitted spectral indices are consider- 
ably lower, over 50% lower for the case with the strongest den- 
sity fluctuations. For this level of turbulence in the background 
plasma the fitted spectral index is always at least a half of the 
source index, indicating the consistent flattening and hardening 
of the spectrum. The only exception to this is for the BW simula- 
tion with the hardest source spectrum (6b = 2.5) where the fitted 
index is 20% higher. This steeper spectrum is because the loss 
of higher energy electrons from the initially harder distribution 
(through the generation of Langmuir waves) dominates over the 
re-acceleration through the plasma inhomogeneities. 



i i 

1.4 - 




2 3 4 5 6 

8 b 

Fig. 7. Power-law index of the electron distribution £tt found 
from fitting f_thick2.pro to the simulated footpoint X-ray 
spectra as a function of the actual initial power-law index 6b- 
The index i5jt is normalised by those found from the fit to the 
beam-only (B, black line) X-ray spectrum. 



The total number of electrons in the source distribution in- 
ferred from the £_thick2 .pro fits is substantially smaller than 
the true values (first panel in Figure |8). The CTTM interpreta- 
tion of these spectra is a considerable underestimate (10 to 1 000 
times) of the number of electrons in the source distribution. It is 
clear in the mean electron spectrum (Figure |4} that including the 
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Fig. 8. Total number of electrons above Eq, 20 keV and 50 keV (left to right) as a function of the power-law index of the initial elec- 
tron distribution 6b, obtained from the fitting of £_thick2 .pro to the simulated footpoint X-ray spectra. The differently coloured 
lines indicate the different simulation setups used. The number of electrons are normalised by those from the fit to the beam-only 
(B, black line) X-ray spectrum. 



plasma inhomogeneity accelerates electrons to higher energies, 
which increases the population at energies well above 20 keV. 
Therefore we also calculated the number of electrons that the 
f _thick2 . pro fit suggests is above 20 and 50 keV in the source 
distribution, as shown in the middle and last panels of Figure 
[8] At best, with steep spectra and the strongest fluctuations con- 
sidered here, the wave scattering can produce a similar number 
of electrons to those found in the CTTM case for the number 
of electrons above 20 keV. For the highest energy electrons the 
situation is considerably better with several simulation setups 
producing number of electrons N{> 50keV) similar to or larger 
than the B case. For the steepest initial spectrum and strongest 
fluctuations the f_thick2 .pro fit overestimates the number of 
electrons in the source by over an order-of-magnitude, about 
a factor of 20. However, it is only a limited set of conditions 
(Zr 1 > 10~ 7 cm _I and 6b > 3.5) that produces more high-energy 
electrons than f_thick2 .pro. 

4. Discussion and conclusions 

Including Langmuir waves driven by the propagating electron 
beam causes major changes in the energy of the electrons and 
produces substantially different X-ray spectra. With no, or low 
levels ZT 1 < 10 _8 cm _I of the density fluctuations in the back- 
ground plasma, the dominant effect is wave-particle interactions 
that decelerate the electrons, which produces a flatter spectrum 
and weaker X-ray emission. If these simulated spectra were as- 
sumed to be caused by the CTTM the number of electrons in 
the source distribution would be substantially underestimated . 
With strong inhomogeneities (L > 10" 7 cm _1 ) in the back- 
ground plasma there is more re-acceleration of the electrons to 
higher energies, resulting in harder (flatter, smaller 6) spectra. 
Interpreting these simulations with the CTTM produces either a 
similar amount or an overestimate (we found up to x20) to the 
number of electrons in the source distribution. Langmuir waves, 
if generated in solar flares, can produce substantial changes in 
the flare-accelerated electron distribution. These effects need to 
be included for a more reliable interpretation of flare HXR spec- 
tra. 

We found for L 1 > lO ^ctrT 1 and 6b > 3.5, the electron 
re-acceleration becomes sharply pronounced, while the low den- 
sity gra dients are insufficient to shift the energy quickly enough. 
Indeed. iRatcliffe et al.l (120121) showed that the strongest acceler- 
ation is achieved when the relaxation time is close to the time 
scale due to the density inhomogeneity. Because we have a fixed 



upper limit to the simulation grid v = 1 15vt the re-acceleration 
of electrons in the hardest source distributions (6 < 3.5) might 
be lost. We are developing a full relativistic treatment of Eqs. (Q]) 
and (0 to study whether the plasma inhomogeneities can have a 
greater effect for this flatter spectral domain. 

Our simulations lack wave-wave interactions. The interac- 
tion of Langmuir waves with ion-sound wav es has been shown t o 
produce additional electron re-acceleration (iKontar et alj|2012h . 
These simulations had no spatial dependence and work is un- 
der way to investigate their role in the ID simulations presented 
here. The interactio n of Langmuir waves with whistler or kinetic 
Alfen waves (e.g. lBian et al.l2 010) might also produce consider- 
able changes to the electron distribution in flares. The density 
fluctuations can effectively change the direction of Langmuir 
waves, and hence depart from the ID model , which could limit 
the ap plication of our simulations. Recently Karlick v & Kontarl 
d2012l) have performed a number of 3D particle-in-cell PIC sim- 
ulations with initially mono-energetic beams and have shown 
that during 3D relaxation a population of electrons appear that 
has velocities exceeding those of the injected electrons. While 
PIC simulations cannot predict the long-term evolution of these 
processes as considered here, the number of accelerated elec- 
trons at the stage of plateau formation closely matches the num- 
bers estimated using ID quasilinear equations. 
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